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Abstract In a second-order r-mode theory, Sa & Tome found that the r-mode oscillation 
in neutron stars (NSs) could induce stellar differential rotation, which leads to a saturation 
state of the oscillation spontaneously. Based on a consideration of the coupling of the r- 
modes and the stellar spin and thermal evolutions, we carefully investigate the influences 
of the r-mode-induced differential rotation on the long-term evolutions of isolated NSs 
and NSs in low-mass X-ray binaries, where the viscous damping of the r-modes and its 
resultant effects are taken into account. The numerical results show that, for both kinds 
of NSs, the differential rotation can prolong the duration of the r-mode saturation state 
significantly. As a result, the stars can keep nearly constant temperature and angular ve- 
locity over a thousand years. Moreover, due to the long-term steady rotation of the stars, 
persistent quasi-monochromatic gravitational wave radiation could be expected, which in- 
creases the detectibility of gravitational waves from both nascent and accreting old NSs. 
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1 INTRODUCTION 

/?-modes in a perfect fluid star with arbitrary rotation arise due to the action of the Coriolis force with 
positive feedback (Andersson 1998; Friedman & Morsink 1998), succumbing to gravitational radiation- 
driven Chandrasekhar-Friedman-Schutz instability (Chandrasekhar 1970; Friedman & Schutz 1978). 
In contrast, the growth of the modes can be suppressed by the viscosity of the stellar matter. Thus 
the r-mode evolution is determined by the competition between the viscous damping effect and the 
destabilizing effect due to gravitational radiation. Based on the conservation of angular momentum, 
a phenomenological model describing the r-mode evolution was proposed by Owen et al. (1998) and 
improved by Ho & Lai (2000). However, since nonlinear effects are ignored in this original version of 
the model, an unbounded growth could lead the modes to an unphysical regime. By putting a saturation 
value for the r-mode amplitude into the model by hand, some authors (e.g., Owen et al. 1998; Levin 
1999; Ho & Lai 2000; Watts & Andersson 2002; Heyl 2002) studied the spin and thermal evolutions 
and gravitational wave radiation of neutron stars (NSs). 

To understand r-modes more deeply and judge their astrophysical implications, it is necessary to 
take into account some nonlinear effects that could give a saturation r-mode amplitude spontaneously 
(e.g., Schenk et al. 2002; Arras et al. 2003; Brink et al. 2004a, 2004b, 2005). As an important nonlinear 
effect, differential rotation induced by r-modes was first studied by Rezzolla et al. (2000, 2001) analyt- 
ically using linearized fluid equations by expanding the velocity of a fluid element located at a certain 
point in powers of the mode amplitude, averaging over a gyration, and retaining only the lowest-order 
nonvanishing term. Soon afterwards, some numerical studies (Stergioulas & Font 2001; Lindblom et al. 
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2001) confirmed the existence of such drifts. More exactly, Sa (2004) solved the fluid equations within 
nonlinear theory up to the second order in the mode amplitude and described the differential rotation 
analytically. By extending the r-mode evolution model of Owen et al. (1998) to this nonlinear case, Sa & 
Tome (2005, 2006) obtained a saturation amplitude of r-modes self-consistently. They also studied the 
early part (millions of seconds after the birth) of the spin evolution of nascent NSs under the influence 
of the differential rotation, but their calculation was not long enough to cover the phase during which 
the effect of the viscous damping of the r-mode operates. In this paper, we would find that the long-term 
spin and thermal evolution of isolated NSs and NSs in low-mass X-ray binaries (LMXBs) can also be 
remarkably influenced by the differential rotation by prolonging the duration of the r-modes. Moreover, 
in view of the prolonged r-modes, it can be accepted that gravitational waves would be continuously 
emitted from both young and accreting old NSs for a long time. 

In the next section, we review the second-order r-mode theory of Sa (2004) briefly. Then, we exhibit 
the coupling thermal, r-mode, and spin evolution equations in Section 3, where some typical numerical 
solutions are given for both isolated and accreting NSs. In Section 4, we estimate the detectability of 
gravitational waves from NSs. Finally, a summary and discussion are given in Section 5. 

2 THE SECOND-ORDER fl-MODES 

For a rotating barotropic Newtonian star, the r-mode solutions of perturbed fluid equations can be found 
in spherical coordinates (r, 6, <fi) at first order in a as (Lindblom et al. 1998), 

S^v r = 0, (1) 

5 {1) v e = a nCil(^j 1 sin^ 1 6»sin(/(/) + cjt), (2) 



t^V = aVLCil[^j sin'" 2 9 cosO cos(Z</> + cjt), (3) 

and at second order in a as (Sa 2004) 

S^v r = 6^v e = 0, (4) 

+a 2 fLAr N - 1 sin"- 1 6, (5) 

where a represents the amplitude of the oscillation, R and O are the radius and angular velocity of the 
unperturbed star, u = -Q(Z + 2)(Z - + 1), C x = (21 - 1)!! V /(2Z + 1)/[2tt(2Z)!Z(Z + 1)], A and 
N are two constants determined by the initial condition. For simplicity, Sa & Tome (2005) suggested 
N = 21 — 1 and redefined A by introducing a new free parameter K as A = ^KCfl 2 (l + 1)R 2 ~ 2L . 
For the most unstable Z = 2 r-mode of primary interest to us, the second-order solution S^v^ shows a 
differential rotation of the star induced by the r-mode oscillation, i.e., large scale drifts of fluid elements 
along stellar latitudes. Using 5^v l and S^v\ the corresponding Lagrangian displacements and 
£( 2 ) 1 can be derived and then the physical angular momentum of the Z = 2 r-mode can be calculated up 
to the second order in a as (Sa 2004; Sa & Tome 2005) 

J r = J (1) + = (4K + 5) a 2 JMR 2 n, (6) 

where J = 1.635 x 10~ 2 and 

= -j P d,e i)i (d t e } + «*v fc ^) dv, (?) 



v(v4 2) +v fe d 2) 



dV. (8) 
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Meanwhile, following Owen et al. (1998) and Sa (2004), we further express the energy of the I = 2 
r-mode by 

E r = - 1 J«0 = (AK + 9) a 2 jMR 2 tf. (9) 

When K = —2, j' 2 ' vanishes and the expressions of J r and E r return to their canonical forms (Owen 
et al. 1998), in other words, the differential rotation disappears. Both the physical angular momentum 
and energy of r-modes are increased by gravitational radiation back reaction and decreased by viscous 
damping, which yields 

A T 9 T 9 T 

(10) 
(ID 

where r g = 3.26^" 6 s, t sv = 2.52 x 10 8 T 9 2 s, and r bv = 6.99 x 10 8 T g - 6 fh 2 s are the timescales of 
the gravitational radiation, shear viscous damping, and bulk viscous damping (for I = 2), respectively 
(Owen et al. 1998), and r v = (r" 1 + t^ 1 ) -1 . Hereafter, the convention Q x = Q/10 X and Cl = 
$1/ \fitGp are adopted in cgs units. These timescales are obtained with a polytropic equation of state as 
p = kp 2 for NSs, with k chosen so that the mass and radius of the star are M = 1AM@ and R = 12.53 
km. The competition between the gravitational destabilizing effect that is dependent on £1 and the In- 
dependent viscous damping effect determines an instability window in the T — CI plane, where a small 
perturbation would grow exponentially due to (t^ 1 — t^ 1 )^ 1 > 0. 

3 EVOLUTIONS OF NSS 

3.1 Thermal evolution equation 

Considering the temperature dependence of the viscosities, we would like to show the thermal evolution 
equation of a NS first before calculating the r mode evolution, which reads (Shapiro & Teuklosky 1983; 
Yakovlev et al. 1999; Yakovlev & Pethick 2004) 



^ = -^(L v +L 7 -H sv ), 



where C v ~ 10 39 T g erg K _1 is the heat capacity of the star. On one hand, the NS can be cooled 
by neutrino and photon energy release, whose luminosities are estimated to be L v w 10 40 Tg erg s" 1 
(for modified URCA process) and L 1 = 4:nR 2 aT^ w 10 35 Xf- 2 erg s _1 , respectively. For the black- 
body luminosity L 7 , the relationship, T s w 3.34 x 10 6 Tg 55 , between the interior (T) and surface 
(T s ) temperatures is used (Gudmundsson et al. 1983). Specifically, the temperature dependence of the 
luminosities indicates that the cooling of the NS at high (> 10 8 K) and low (< 10 8 K) temperatures would 
be dominated by neutrino and photon emissions, respectively. On the other hand, the shear viscous 
dissipation of r-modes can convert a part of the oscillation energy into heat energy gradually. Using the 
shear viscous damping timescale, the rate of this energy conversion can be estimated by 

IE 

H sv = — -= 2.0 x 10 43 (4X + 9)a 2 Tg 2 Cl 2 ergs" 1 . (13) 

Tsv 

For the very early ages of a nascent NS, during which this heating effect is much weaker than the 
neutrino cooling effect yet, an approximative temperature evolution can be solved from Eq. ( fT2l as 
T = Tifi+t/tc) -1 ' 6 , where Ti is the initial temperature and t c « (20/Xf 10 )s. However, as the r-modes 
increase, the cooling of the star would be resisted effectively by the heating effect, as demonstrated by 
some previous studies (e.g., Zheng et al. 2006). 
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3.2 Isolated NSs 

A simple phenomenological model for r-mode evolution was proposed by Owen et al. (1998) first and 
further improved by Ho & Lai (2000) based on a consideration of angular momentum conservation. For 
a normal NS with a strong magnetic field (~ 10 10 ~ 12 G), besides the braking effect due to gravitational 
radiation, the spindown of the star resulting from magnetic dipole radiation should also be taken into 
account. So, we ought to write the decrease of the total angular momentum of the star as (Owen et al. 
1998; Ho & Lai 2000; Sa & Tome 2005) 

dJ _ 3a 2 JMR 2 Q m 

(It Tg T m 

where r m = 1.35 x 10 9 B^ 2 (Q/\/t:Gp)^ 2 s is the magnetic braking timescale and / = IMR 2 with I = 
0.261 is the moment of inertial of the star. Due to the /--mode oscillation, the total angular momentum 
of the star could be separated into two parts, i.e., J = IVl + J r . Then, Eqs. ( [Tol l and ( fT4l yield 



da 
~dt 



\ + \{K + 2)Qa 7 



•i 



\ + \{AK + b)Qa 7 



-+2^ (15) 



dn 8 , o ft 2 , , 9 ft O 

— = --(K + 2)Qa 2 - + -(4K + 5)Qa 2 , (16) 

at 6 Tg 6 T V 

Tm 

where Q = 3J/2I = 0.094. During the very early ages of nascent NSs when T g <C (t u , r m ), the viscous 
and magnetic terms in the above two equations can be omitted. Combining this simplification with the 
analytical temperature T = Tj(l + i/i c ) -1 / 6 , Sa & Tome (2005, 2006) obtained an analytical solution 
of Eqs. ( Tl3T > and ( TTST ) for t < O.lyr. For convenience, their analytical solution can also be expressed by 
two asymptotic functions as follows (Sa & Tome 2006): 

( an exp {t/T g ,i) , for<<i Q 

m ~ \0.63(t/r ff . 4 )- 1/5 , for t > t a (W 

where a,i and are the initial r-mode amplitude and angular velocity, respectively, and r Si , = 
3.26f2^~ s. The transition time t a f» [521 — 18.51n(_RT + 2)]s corresponding to the amplitude of 
a{t a ) = [12(K + 2)Q]" 1 / 2 is determined by the condition d 2 a/dt 2 = (Sa & Tome 2006). 

However, as the temperature and angular velocity decrease, the viscous damping timescale would 
become to be comparable to the gravitational radiation timescale. Therefore, it is necessary to com- 
pletely solve the coupling Eqs. ( fT2b . (1151) . and (fTST i in order to depict the long-term history of NSs. For 
different values of K (> —5/4), we show some numerical evolution curves of the r-mode amplitude in 
Figure 1 . As indicated by the thin solid lines, the two increasing segments of the evolution curves can be 
fitted by Eq. ( fT71 ) well, i.e., the amplitude increases rapidly first and then gradually reaches a saturation 
value. About one tenth year later after the birth of the stars, the growth of the r-mode would be stopped 
and instead, the amplitude nearly keeps constant until an extremely fast decay due to (r ~ 1 — t~ 1 ) ~ 1 < 0. 
The higher the value of K, the longer the duration of this plateau phase. 

In order to exhibit the influence of the differential rotation on the r-mode evolution, for an example, 
we plot the r-mode evolution curves for K =100 (differential rotation case) and —2 (non-differential 
rotation case) in Figure 2(a) for a comparison. As mentioned above, the non-differential rotation model 
(K = 2) is incapable of determining a saturation amplitude. So, in the case of K = —2, we put an 
effective saturation amplitude by hand, which is taken to equal the one calculated from the contrastive 
differential rotation case for consistency (e.g., a sat = 1.1 for K = 100). Correspondingly, Figures 
2(b) and 2(c) show the temporal evolution of the stellar angular velocity and temperature, respectively, 
for both K = 100 and —2. Especially, for the differential rotation case, we divide the stellar evolution 
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Table 1 Different phases of the evolution of a young NS during the r-mode oscillation. The 
temporal behaviors of a, CI, and T are listed for every phase. The coefficient a = ^(K + 

2)Qal 



Phases 


I 


II 


Ilia 


Illb 


IV 


V 


a(t) oc 


exp(t/r 9 ,i) 


exp(t/r g ,i) 


t 1/w 


t 1/w 


i° 


decrease 


fi(t) oc 


t° 


1 - aexp(2t/r 9ii ) 


t -i/s 


t -l/5 


t° 


increase 


T(t) oc 


t -l/6 


t -l/6 


t -l/B 


t° 


t° 


decrease 



during the r-mode oscillation into six phases (denoted by I- V) roughly, the temporal behaviors of which 
are listed in Table 1 . Within phase IV, the slow changes of fl and T make the timescales r g and t v vary 
slowly. Thus, the simultaneous r-mode oscillation can maintain steady for a long period. Comparing the 
differential with non-differential rotation cases, we can find that: (1) The differential rotation obviously 
strengthens the gravitational braking effect for t < 0.1 yrs (phases II and III). However, subsequently, 
from one tenth to a few thousand years (phases IV and V), the spindown of the star due to gravitational 
radiation would be held back effectively by an angular momentum transfer from J r to IQ,, although dur- 
ing this time the r-mode always stays in the saturation state. Due to the existence of this angular velocity 
plateau (i.e., dVl/dt ~ 0; phase IV), the star is expected to emit a quasi-monochromatic gravitational 
wave persistently. (2) The obvious difference in the temperature plateaus between the cooling curves 
with K = 100 and —2 indicates that the heating effect due to r-mode dissipation is also strengthened 
dramatically by the differential rotation. As a result, the NSs with differential rotation can keep a high 
constant temperature for a few thousand years. In view of the nearly constant temperature and angular 
velocity within phase IV, it is easy to understand the appearance of the steady r-mode saturation state. 
Finally, we also show the evolution trajectories of an isolated NS for K = 100 and —2 in the T — O 
plane in Figure 3, where the six evolution phases defined for the differential rotation case are labelled. 
Especially, within phase V, a self-spinup of the differential-rotation star can be seen clearly. In addition, 
phase IV is marked by a solid circle, where a quasi-monochromatic gravitational wave could be emitted 
for a few hundred years (see Sect. 4). 

To summarize, during the early part of the r-mode evolution (phases I, II, and III), the rotation 
energy of the star (|IQ 2 ) is converted into the oscillation energy, the internal energy, and the energy 
of gravitational waves. In contrast, during the late part (phases IV and V), the energy deposited in the 
r-modes would be released gradually via heating the star and accelerating the stellar rotation due to 
viscosity. Moreover, this spin-up effect could be stronger than the gravitational braking effect. 



3.3 NSs in LMXBs 

For NSs in LMXBs, whose magnetic fields are usually found to be relatively weak (~ 10 8-9 G), their 
angular velocity could be increased by accreting materials from their companian star. Then, the evolution 
of the stellar angular momentum would be determined by the competition between the gravitational 
radiation and accretion as (Levin 1999; Zhang & Dai 2008) 

g = _3a»JM#n +J . 

(it Tg 

where M is the accretion rate and the velocity of the accretion disk is assumed to equal the Keplerian 
velocity f2#. Combining Eqs. ( flUt and ( fT9l , we can get 
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where r a = M/M is defined as an accretion timescale. 

We plot in Figure 4 the evolution trajectories of an accreting NS in the T — SI plane for K = 100 
and —2. Different from the case of the isolated NS shown in Figure 3, the accreting star can be spun 
up by accretion significantly rather than spun down by magnetic dipole radiation in old age (~ 10 5-6 
yrs). Especially, if the accretion rate is high enough, cyclic evolution could be found (black lines). 
This is qualitatively consistent with the results in Levin (1999) and Heyl (2002). However, for M = 
10 _8 M Q yr" 1 specifically, we do not obtain the cycle but Levin (1999) did. There are two reasons for 
this difference: (1) In the calculations of Levin (1999), an effective shear viscous damping timescale 
t S v = 1-03 x 10 6 Tg s was taken by hand in order to fit the observed data, whereas we adopt a theoretical 
value of t sv = 2.52 x 10 8 T g 2 s from Owen et al. (1998); (2) The cooling effect due to the thermal 
radiation, which can effectively pull the star away from the r-mode instability window in the T — SI 
plane, was ignored in Levin (1999). 

The temporal behaviors of a, SI, and T within one cycle are exhibited in Figure 5. In order to show 
the detailed features of the cycle clearly, the time-axes in the left- and right-hand panels of Figure 5 
are drawn on normal and logarithmical scales, respectively. To be specific, the left-hand panel shows 
that the period of the cyclic evolution is shortened by the differential rotation mildly (4.5 x 10 5 yrs 
vs 5.6 x 10 5 yrs), and the right-hand panel indicates that the duration of the r-mode oscillation within 
one cycle is prolonged significantly (3,900 yrs vs 65 yrs). Similar to the early evolution of young NSs 
shown in Figure 2, the evolution during the r-mode oscillation within one cycle of the accreting NSs can 
be divided into five phases. A comparison between Figures 2 and 5 shows that the temporal behaviors 
of phases Illb, IV, and V of young and accreting old NSs seems to be nearly identical except for their 
durations. This indicates that the isolated young and accreting old NSs may be able to produce some 
same astrophysical phenomena, e.g., self-spinup and persistent quasi-monochromatic gravitational wave 
radiation. 



4 DETECTABILITY OF GRAVITATIONAL WAVES FROM THE /WVIODE 

Using the obtained r-mode amplitude and angular velocity, we can estimate the amplitude of the emitted 
gravitational waves as follows (Owen et al. 1998; Sa & Tome 2006): 



\h(t)\ = 1.3 x lQ- 24 a{t) 



~Sl(t)~ 


3 


^20Mpc\ 


[n K _ 




V d L ) 



(22) 



where is the distance of the star. Then, the frequency-domain gravitational wave amplitude [i.e., the 
Fourier transform of h(t), h(f) = e 2mft h(t)dt] can be calculated by (Owen et al. 1998; Sa & 
Tome 2006) 

where / = 2Sl/(3ir) is the frequency of the gravitational waves. In order to analyze the possibility of 
detecting the gravitational waves with laser interferometer detectors LIGO and Virgo, in Figure 6 we 
compare the characteristic amplitude of the signal, h c (f) = f\h(f)\, with the rms stain noise in the 
detectors, h Tms (f) = y fSh(f), for both isolated (left-hand panel) and accreting (right-hand panel) 
NSs. For the noise spectral density of the detectors, Sh(f), some approximative expressions can be 
found for LIGO, Virgo, and advanced LIGO in Sa & Tome (2006). 

On one hand, as found by Sa & Tome (2006), the spike of h c (f) at / max = 2SIk /(3tt) that was 
predicted by Owen et al. (1998; see the thick dashed lines in Figure 6) disappears under the influence of 
the differential rotation, and the numerical results of h c (f) for / > 100 Hz in Figure 6 can be fitted by 
the following analytical expression perfectly 
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Table 2 Signal-to-noise ratio of gravitational wave detections for different detectors, different 
values of K, and different detection duration for an isolated NS at cLl = 20 Mpc. The begin 
of the detection is set at the birth of the star. 





L1GO 


Virgo 


advanced L1GO 


tdct — to 


K = 1 


10 


100 


1 


10 


100 


1 


10 


100 


0.3 yrs 


0.65 


0.33 


0.11 


0.51 


0.26 


0.01 


9.18 


4.62 


1.59 


1 yrs 


0.67 


0.34 


0.12 


0.54 


0.27 


0.01 


9.63 


4.90 


1.69 


10 yrs 


0.81 


0.44 


0.15 


0.77 


0.42 


0.15 


13.53 


7.46 


2.63 


30 yrs 


0.91 


0.55 


0.20 


1.06 


0.62 


0.22 


18.60 


10.98 


3.93 



On the other hand, surprisingly, a new remarkable spike emerges within the range of ~ 60 — 90 Hz, 
where the approximative analysis in Sa & Tome (2006) is inapplicable. From Figures 2 and 5 we know 
that, during phase IV, the angular velocity of the star can nearly keep constant (i.e., \df/dt\ — > 0) 
over a few hundred years, while the r-mode stays in the saturation state all the time. As a result, a quasi- 
monochromatic gravitational wave could be emitted from both young and accreting old NSs, which lasts 
a few hundred years. Additionally, for accreting NSs, another weaker spike at ~ 220 Hz is predicted 
due to the existence of Phase I'. 

Using matched filtering, the power signal-to-noise ratio (S/N) 2 of a detection from to to tdet is 
given by (Owen et al. 1998; Sa & Tome 2006) 

w 4 k„ s (/(d ' ( } 

where to is determined by the begin of the detection. In Table 2, we list some values of S/N with 
different tdct and K for an isolated NS by setting to at the birth of the star. Since the spike of h c (f) 
within ~ 60 — 90 Hz would appear about 0.3 yrs later after the rising of the r-modes, the signal-to-noise 
ratio obtained from a long-term detection could be much higher than that from a short-period detection 
(i.e., the case focused in Sa & Tome 2006). 

5 SUMMARY AND DISCUSSION 

A second-order r-mode theory was developed by Sa & Tome (2004; 2005). This theory predicts that 
r-mode oscillation could induce differential rotation in neutron stars, which can determine a saturation 
amplitude of the r-mode spontaneously. In the framework of this theory, we investigate the long-term 
spin and thermal evolutions of isolated NSs and NSs in LMXBs. In our calculations, the effects of heat- 
ing due to r-mode dissipation, gravitational and magnetic braking, and accretion are taken into account. 
Our results show that, to a certain extent, the linear r-mode evolution model using an artificial saturation 
amplitude can describe the basic features of the evolution of NSs qualitatively, but predicts an obviously 
underestimated r-mode duration. By considering the differential rotation, we may obtain a slight self- 
spinup and an enhanced temperature plateau for NSs. Especially, due to the effective angular momentum 
transfer from J r to lil, the spindown of the NSs can be stopped for a few hundred years, whereas the 
gravitational radiation still exists during this period. As a result, long-lasting quasi-monochromatic grav- 
itational wave radiation is predicted, which increases the detectability of gravitational waves from both 
nascent and accreting old (~ 10 5 ~ 6 yrs) NSs. 

In this paper, we adopt a very simple NS model just in order to find the influences of nonlinear 
effects on the evolution of NSs qualitatively. However, generally speaking, NSs are probably hybrid 
stars or even strange quark stars. The former undergoes a deconfinement transition from neutron matter 
to quark or hyperon matter (Glendenning 1997; Pan et al. 2006), and the latter consists of nearly pure 
quark matter (Alcock et al. 1986; Zheng et al. 2006). It is a demanding task to study the r-mode evolution 
in a more realistic NS model. Especially, for a hybrid star that contains a quark or hyperon core, since 
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the direct Urea is triggered and superfluidity exists, the star without reheating would have a very low 
temperature that is inconsistent with the observational data. So, it may be helpful for enhancing the 
temperature of the hybrid stars to consider the heating effect due to r-mode dissipation. 
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Long-term evolution of neutron stars 




Fig. 1 Evolution of r-mode amplitude of an isolated NS with a magnetic field B = 10 G for 
different values of K (thick lines). The thin solid lines are given by the asymptotic functions 
shown in Eq. ( fTTI i. The initial values of the r-mode amplitude, angular velocity, and tempera- 
ture are taken to be = 10~ 6 , fij = = and T, = 10 10 K, respectively, where 
VLk is the Keplerian angular velocity at which the star starts shedding mass at the equator. 
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Fig. 2 Evolution curves of a, f2, and T of an isolated NS with a magnetic field £> = 10 12 G for 
K = 100 (solid lines; differential rotation case) and K = —2 (dashed lines; non-differential 
rotation case). The initial conditions are the same to that in Figure 1. The evolution during the 
r-mode oscillation is divided into several phases (denoted by I-V) by the vertical dotted lines. 
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Fig. 3 Evolution trajectories of an isolated NS with B = 10 12 G in the T — Q, plane for 
K = 100 (solid line; differential rotation case) and K = —2 (dashed line; non-differential 
rotation case). The initial conditions and the meaning of phases I-V are the same to that in 
Figure 1 . The shaded region exhibits the r-mode instability window. 
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Fig. 4 The same to Figure 3 but for a NS with B = 10 8 G in a LMXB. The black and grey 
lines correspond to the accretion rates M = lO^M^yr -1 and lO _8 M0yr _1 , respectively. 
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Fig. 5 Evolution curves of a, 0, and T during the cyclic evolution of a NS with B = 10 8 G 
in a LMXB for K = 100 (solid lines; differential rotation case) and K — —2 (dashed lines; 
non-differential rotation case). The begin of the cycle is set at point A that is marked in Figure 
4, and the age of the star at point A is denoted by t C ycic,o- 
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Fig. 6 A comparison of the characteristic amplitude of gravitational waves for different values 
of K (thick lines) with the rms strain noise in the detectors (thin solid lines). The details of 
the spike of h c within ~ 60 — 90 Hz is shown in the insert panel. 



